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Abstract. Recent work giving a classification of kinematic and vorticity conservation laws 
of compressible fluid flow for barotropic equations of state (where pressure is a function only 
of the fluid density) in n > 1 spatial dimensions is extended to general non-isentropic equa- 
tions of state in which the pressure is also a function of the dynamical entropy (per unit 
mass) of the fluid. Two main results are obtained. First, we find that apart from the fa- 
miliar conserved integrals for mass, momentum, energy, angular momentum and Galilean 
momentum, and volumetric entropy, additional kinematic conserved integrals arise only for 
non-isentropic equations of state given by a generalized form of the well-known polytropic 
equation of state with dimension-dependent exponent 7 = \ -\- 2/n, such that the propor- 
tionality coefficient is an arbitrary function of the entropy (per unit mass). Second, we show 
that the only vorticity conserved integrals consist of a circulatory entropy (which vanishes 
precisely when the fiuid flow is irrotational) in all even dimensions. In particular, the vor- 
ticity integrals for helicity in odd dimensions and enstrophy in even dimensions are found 
to be no longer conserved for any non-isentropic equation of state. 



1. Introduction and summary 

The mathematical study of n-dimensional fluid flow has attracted rising interest in the 
past few decades, encompassing work on symmetries and conservation laws j211 [131 E]; 
Hamiltonian structures [25j, Casimir invariants |22l [121 [15] , ^"^^ other group-theoretic aspects 
of the n-dimensional Eulerian fluid equations [3 [El [9] . Further results in the special case of 
n = 1 dimension appear in [231 1211 [20] • 

In a recent contribution [5], we undertook a systematic study of local conservation laws for 
the Euler equations governing isentropic compressible fluid flow in n > 1 spatial dimensions, 
where the pressure of the fluid is a function of the fluid density as given by some barotropic 
equation of state, while the entropy (per unit mass) of the fluid is constant throughout 
the fluid domain. Our results fully settled the problem of finding all n-dimensional local 
continuity equations in two cases of primary interest: kinematic conservation laws, like 
mass, momentum and energy, for which the conserved density and spatial flux involve only 
the fluid velocity, density and pressure, in addition to the time and space coordinates; and 
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vorticity conservation laws, such as helicity in three dimensions as well as circulation and 
enstrophy in two dimensions, where the conserved density and spatial flux have an essential 
dependence on the curl of the fluid velocity. These two classes of conservation laws comprise 
all of the local continuity equations known to-date for isentropic compressible fluid flow in 
n > 1 dimensions. 

The present paper extends this work to non-isentropic adiabatic compressible fluid flow 
in which the entropy (per unit mass) is conserved only along streamlines and the pressure is 
given by an equation of state in terms of both the fluid density and entropy. In particular, we 
explicitly derive all n-dimensional kinematic and vorticity conservation laws, including any 
that are admitted only for special non-isentropic equations of state or in special dimensions 
n > 1, with the conserved densities and spatial fluxes allowed to depend on the entropy (in 
addition to the previous fluid variables). 

In section [21 we begin by reviewing the formulation of necessary and sufficient equa- 
tions for directly determining the conserved densities admitted by the Euler equations for 
n-dimensional non-isentropic compressible fluid flow. We next verify that the physically fa- 
miliar kinematic conserved integrals for mass, momentum, energy, angular momentum and 
Galilean momentum in compressible fluid flow with a general isentropic equation of state 
remain conserved for the n-dimensional non-isentropic Euler equations. By solving the de- 
termining equations for kinematic conserved densities in n > 1 dimensions, we then show 
that the only additional conserved integrals consist of volumetric entropy in a generalized 
form, plus dilational and similarity type energies arising for polytropic equations of state 
where the pressure is proportional to a particular dimension-dependent power 7 = 1 + 2/nof 
the density, with the proportionality coefficient given by an arbitrary function of the entropy 
(per unit mass). 

In section [3l we first consider the odd-dimensional helicity integral and the even- 
dimensional enstrophy integral both of which are known to be conserved in n-dimensional 
isentropic fluid flow [151 E] • By examining how the corresponding local vorticity conservation 
laws depend on the equation of state of the fluid, we show that the helicity and enstrophy 
integrals are no longer conserved for any non-isentropic equations of state. Next we solve 
the determining equations to see if there are any vorticity conserved densities in n > 1 
dimensions admitted by the non-isentropic Euler equations. As a main new result, this clas- 
sification yields a non-trivial conserved integral given by an underlying local conservation 
law describing a circulatory entropy (which vanishes whenever the fluid is irrotational) in all 
even dimensions. 

Next, in section HI we classify kinematic and vorticity constants of motion on n-dimensional 
domains that move with the fluid. Such constants of motion arise from non-trivial local 
conservation laws such that the net flux through the moving boundary of the domain is 
zero. The main result of this classification is to show that, firstly, there are no vorticity 
constants of motion and, secondly, the only kinematic constants of motion are given by 
the mass integral and the volumetric entropy integral defined on any moving domain in the 
fluid. We also investigate possibilities for constants of motion defined by a conserved integral 
on the boundary of a moving domain, in analogy with Kelvin's circulation integral in the 
case of isentropic fluid flow in = 2 dimensions. Such constants of motion correspond to 
conservation laws satisfying the condition of zero net flux on the moving boundary but such 
that the conserved density in the moving domain is locally trivial (i.e. having the form of a 
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spatial divergence). From our classification proofs for kinematic and vorticity conservation 
laws, we find that there are no moving-boundary constants of motion for non-isentropic 
compressible fluid flow in n > 1 dimensions, in contrast to the generalized circulation integral 
which exists in the isentropic case for all even dimensions as found in 

In section El we then use a Hamiltonian formulation of the non-isentropic compressible 
Euler equations to classify all Hamiltonian symmetries corresponding to the kinematic and 
vorticity conservation laws in n > 1 dimensions. 

In section ini based on the well-known relationship between the non-isentropic compressible 
Euler equations and the equations of gas dynamics, we state a complete classification of 
conservation laws of kinematic and vorticity forms for gas dynamics in n > 1 dimensions. In 
this classification, the fluid equation of state gets replaced by a state function that specifies 
both the speed of sound in the gas and the thermodynamic energy density of the gas in terms 
of the gas pressure and density. We thereby obtain conserved integrals for mass, momentum, 
energy, angular momentum and Galilean momentum, plus volumetric entropy and, in even 
dimensions, circulatory entropy, holding for all state functions. The only distinguished state 
function for which the gas dynamics equations admit additional conserved integrals is shown 
to be the polytropic case, where the state function is proportional to the pressure, with a 
particular dimension-dependent proportionality coefficient 7 = 1 + 2/n. 

Finally, in section [TJ we summarize our classification of fluid flow conservation laws in 
index notation and also list the multipliers that correspond to each non-trivial conserved 
density. Proofs of the main classification theorems are given in two appendices. 

An interesting question left open for future work is to extend these classifications to 
conservation laws that depend on derivatives of the curl of the velocity or more generally on 
higher derivatives of the velocity itself. 



2. Non-isentropic compressible Euler Equations 

Compressible non-isentropic fluids in R" are described (in the absence of external forces) by 
a generalization of the Euler equations such that the equation of state for pressure becomes 
a function of the entropy in addition to the density: 

p = P(p,S), (2.1) 

9tu + u ■ Vu + -Vp = 0, (2.2) 
P 

dtp + V- (pu) = 0. (2.3) 

Here u(t,x) is velocity, p{t,iL) is density, and 5'(t,x) is entropy (per unit mass) of the fluid 
which is conserved along streamlines 

dtS + VL-VS = 0. (2.4) 

In particular when S =const. or dP/dS = 0, the system fl2.ip - fl2.3p reduces to the case of 
isentropic fluid flow studied in |5] where 

P = Pip) (2.5) 

so thus the non-isentropic case is characterized by having dP/dS ^ with 5* being a dy- 
namical variable which is non-constant across streamlines. 
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We will study local conservation laws of the Eulerian system f l2.ip -( l27i|) . Similar to the 
isentropic case, conservation laws of non-isentropic fluid flow are described by a local conti- 
nuity equation 

DtT + ■ X = (2.6) 

holding formally for all solutions of fl2.ip - fl2.4p . where the conserved density T and spatial 
flux X are some functions of t, x, u, p, S, and their x- derivatives. Here Dt and denote 
total time and space derivatives respectively. In integral form, the continuity equation (12. 6p 
is equivalently given by 

^ / T^x = - [ X ■ nd"-V (2.7) 
dt Jv Jdv 

where V is any spatial domain in M" through which the fluid is flowing, fi is the outward unit 
normal and d"'~^a denotes the surface element of the domain boundary dV. Alternatively, 
conservation laws (12. 6p and (12. 7p can be formulated by considering a spatial domain V(t) 
that moves with the fluid. Then the spatial flux through the moving boundary dV{t) is 
^ = X — Tu which is related to the conserved density T by the transport equation 

AT + u • D^T = -(V • u)T - A ■ I (2.8) 

where A + u-Dx represents the total convective (material) derivative and V-u represents the 
expansion or contraction of an infinitesimal volume moving with the fluid. The corresponding 
integral form of a fluid conservation law in a moving domain is thereby expressed as 

^ / rc/"x = - / I ■ nrf"- V 
di Jv{t) Jdv{t) 

in terms of the moving- flux ^ through the domain boundary dV{t). 

For determining conserved densities T, we note the local continuity equation (12. 6 p shows 
that VtT must have the form of a spatial divergence -Dx ■ (~X), where Vt is the total time 
derivative evaluated on solutions of (I2.ip -( l2^ . Hence, necessary and sufficient equations 
jlU] for finding conserved densities T are given by 

E^{VtT) = E,{VtT) = EsiVtT) = (2.10) 

where E^, Ep, Es are spatial Euler operators with respect to u, p, S. To-date, a complete 
classification of all conserved densities T(t, x, u, p, S", Vu, Vp, VS", . . .) has not appeared in 
the literature. 

Throughout, on M", V will denote the gradient operator, while ■ and A will respectively 
denote the Euclidean inner product and exterior (antisymmetric) product. 



(2.9) 



2.1. Classification of kinematic conservation laws. For n-dimensional compressible 
fiuid fiow, we first write down the well-known kinematic conservation laws [HI [S] for mass, 
momentum, angular momentum, Galilean momentum, and energy in terms of the fiuid 
velocity u, density p and pressure p, on any spatial domain V{t) C transported in the 
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fluid: 



d 
dt 
d 
dt 
d_ 
dt 
d 
dt 
d 
dt 



pd'^x = 0, (2.11) 

V{t) 

pud^x = - I pud^'-^a, (2.12) 

V{t) JdV{t) 

puAxc/"x= / pxAnrf"-V, (2.13) 

V{t) JdV(t) 

p{tu - x)rf"x = - / tphd'^-^a, (2.14) 

V{t) JdV{t) 

p(||u|2 + e)rf"a; = - / pu ■ nrf"" V. (2.15) 

V{t) JdV(t) 



= p~^p. (2.16) 

D = const. 



Here e is the internal (thermodynamic) energy density defined in terms of p and p by 

de 
dp 

We note that equation fl2.16p is equivalent to the familiar thermodynamic relation pTl [26] 
d{e + p/p) = {1/ p)dp holding for any adiabatic infinitesimal change dS = Q va. the state of 
the fiuid, where e + p/ p is the enthalpy. 

Each of these conservation laws can be readily verified to hold using just the dynamical 
Euler equations fl2.2l) - fl2.3p for u and p, plus the transport equation (12.41) for S in the case of 
the energy conservation law. As a consequence, the form of the conserved mass, momentum, 
angular momentum, and Galilean momentum is independent of the equation of state for p, 
and hence the conserved integrals (12. lip . (I2.12p . (12.130 . (12.140 are valid for non-isentropic 
fiuid fiow. In contrast, the form of the conserved energy (12.150 explicitly depends on the 
equation of state for p in terms of both p and 5* through the relation (12.160 for the internal 
energy density. In particular, for fiuid fiow with a non-isentropic equation of state (12. ip . the 
internal energy density is given by 

e{p,S) = j p-'Pip,S)dp. (2.17) 

Note this defines e only up to an integration constant given by an arbitrary function of 
S. Such a function f{S) contributes a term of the form pf{S) to the energy density E = 
p(||up + e) of the fiuid, and as a result, conservation of energy yields the conservation law 



d_ 
dt 



pf{S)d''x = 0. (2.18) 

V{t) 



In the specific case f{S) = S, this conservation law (I2.18P states that the volumetric entropy 
is conserved in any spatial domain transported by the fiuid, 

d 



dt 



pSd^x = (2.19) 

V{t) 



which is easy to check directly from the dynamical equations (I2.3p - (l2.4p for p and 5*. The 
general case of the conservation law (I2.18P can be understood to arise from conservation of 
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volumetric entropy fl2.19p by the freedom to functionally redefine the entropy S" to S" = f{S), 
preserving the form of the entropy transport equation ( 12.40 . 

The only other kinematic conservation laws known for compressible fluid flow in n > 1 
dimensions are the following generalized energies [HI [S] : 

^ [ {tE- \pvL ■ x)rf"x = - / p(tu - ix) ■ nc/"-V, (2.20) 
dt Jv(t) Jdv(t) 



d 
dt 



/ (t2^-tpu-x+ip|xnd"a; = - / pt(tu - x) ■ nrf"- V, (2.21) 

Jv(t) JdV{t) 

each holding for a polytropic equation of state 

p = K = const. (2.22) 

where 

E = p(i|u|^ + e) = ip|up + inp (2.23) 

is the polytropic energy density. 

We now settle the natural questions of whether these generalized energies (I2.20l) - fl2.2ip 
hold for any non-isentropic equations of state, and whether the non-isentropic compressible 
fluid equations fl2.ip - (l2.4p admit any additional conservation laws of kinematic form 

T(t,x,u,p,5). (2.24) 

Theorem 2.1. (i) For compressible fluid flow with a general non-isentropic equation of state 
Ii2. the admitted kinematic conservation laws 1^2. 24\ ) in any dimension n > 1 comprise a 
linear combination of mass Ii2.11\) . momentum Ii2.12\) . angular momentum Ii2.13\) . Galilean 
momentum 1^2.14^ , energy Ii2.15\) - ^2.1'7\ ), and generalized entropy Ii2.18\) . (ii) Modulo a con- 



stant shift in the pressure, the only equations of state for which additional conservation laws 
^2. 24\ ) arise for non-isentropic compressible fluid flow is the polytropic case 

p = k{S)p^+^/'' (2.25) 

given in terms of an arbitrary function of the entropy, k,{S). The admitted conservation 
laws consist of the similarity energy Ii2.20\) and the dilational energy Ii2.21\) where e(p, S) = 
^K{S)p^f'^ is the internal energy density and E = p(||up + e) = |p|up + |p is the polytropic 
energy density. 

We give the proof of this classification theorem in Appendix El 

3. VORTICITY CONSERVATION LAWS 

We will start by examining conservation of the odd-dimensional helicity integral and the 
even- dimensional enstrophy integral for compressible fluid flow in tt, > 1 spatial dimensions. 
These integrals involve the spatial orientation tensor e contracted into products of the curl 
of the fluid velocity [5] 

a; = VAu. (3.1) 

The dynamical equation for uj is given by the curl of the Euler equation (12. 2p for u, combined 
with the identity u ■ = u ■ Vu — |V(u ■ u), which yields 

dtU) = V A{u:-u) + (T, V A w = 0, (3.2) 
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where 

a = Vp"^ A Vp, V A cr = 0. (3.3) 

We note that in the case of an isentropic equation of state (12.51) . this antisymmetric tensor 
(13. 3p vanishes identically, since Vp = P'{p)Vp is proportional to Vp. 
In odd dimensions n = 2m + 1 > 3, the integral 

/ u ■ -ccrrf'^x (3.4) 

Jv{t) 

defines the fluid helicity on any spatial domain V{t) C M" transported in the fluid, where 

z:^ = e ■ ((V A u) A • • ■ A (V A u)) = *(a;'") (3.5) 

" V ' 

(n — l)/2 times 

denotes the vorticity vector of the fluid. (Here * is the Hodge dual operator applied to 
the rank n — 1 = 2m skew-symmetric tensor cj™.) Since u: is curl-free, and the tensor e 
is constant and skew-symmetric, the vorticity vector (13.51) is divergence-free and obeys the 
transport equation 

dt^ = -V ■ {x^ Au)+m(T -W, W = *(a;'"-^), (3.6) 

with 

V ■ = *(V A a;™) = m(V A a;) ■ ^iuj""-^) = (3.7) 

obtained from (13. 2p and (13. 5p . Now for evaluating the time derivative of the helicity integral 
(13. 4p . we first use equations (13. 6p and (12. 2p to derive 

dt{u ■ Tu) = V ■ (-u(u ■ Tu) + ^-ccrlul^) - p-^vj ■ Vp - mW ■ (u A cr). (3.8) 

We next rewrite the pressure gradient term in (13. 8p by means of the identities W ■ ( V A u) = 
W • a; = w and V ■ W = (m - 1)(V A a;) ■ *{u}"'-^) = 0, yielding 

p-^-cc7 ■ Vp = V ■ (W ■ (u A Vp)p'^) + W ■ (u A <t) (3.9) 

via integration by parts. We thus obtain 

dtiu • -CC7) + V ■ (u(u ■ -CC7) - itx7|u|2 + W ■ (u A Vp)p-^) = -(1 + m)W ■ (u A <t) (3.10) 

which has the form of a local continuity equation (12.61) up to the term proportional to cr on the 
right hand side. This verifies conservation of helicity in the case of isentropic compressible 
fluid flow. For the general case of non-isentropic compressible fluid flow, helicity will be 
conserved if and only if, for a given equation of state, the non-vanishing cr term on the right 
hand side in equation (I3.10p reduces to a total spatial divergence. Necessary and sufficient 
conditions for this to occur are that the spatial Euler operators Ep, Eg, E^ annihilate the cr 
term. 

To proceed, for a general non-isentropic equation of state p = P{p,S), we note Vp = 
PpVp + PsVS implies cr = -p^'^PgVp AVS = -Ve^ A VS with ^ 0, where e is the 
internal (thermodynamic) energy density (I2.17P defined in terms of P{p,S). Hence, the 
antisymmetric tensor (13.31) in the non-isentropic case is given by the curl 

cr = -V A (esV^). (3.11) 

This allows us to rewrite the right hand side in equation (I3.10p 

W ■ (u A cr) = V ■ (W ■ (u A VS)es) - e^V ■ (v^S) (3.12) 
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via integration by parts. Similarly, we can rewrite the divergence term 

V ■ (W ■ (u A Vp)p-^) = V ■ {vu{p'^p + e) - esW ■ (u A VS)) (3.13) 

by means of the identity 

Ve = epVp + esVS = -pVp'^ + esVS = esVS + p-^Vp - V{p-^p) (3.14) 

which follows from the thermodynamic relation (12.161) . The conservation equation (13.101) 
thereby reduces to 

DtT + ■ (uT) + • ^ = (1 + m)esV ■ (w5) (3.15) 

with 

T = u--cc7, ^ = {p-^p + e-l\u\'^)vj + mesW ■ {VS Au). (3.16) 

We now observe the right hand side in equation (I3.15P fails to be a spatial divergence when 
the fluid flow is non-isentropic, since Pg ^ implies Ep^esV ■ (zuS) = CpsV ■ {'cuS) = 
p~'^Ps'ccr ■ VS ^0. As a result, the helicity integral (13. 4p obeys 

^[ u-zj7rf"x= / (i|u|2 -p-ip-e)-cc7- nd"-V + m / egW ■ (V5 A u A n)rf"- V 
Jv(t) Jav(t) Jdv(t) 

+ (l+m) ! esV ■ {'ujS)<rx (3.17) 
Jv(t) 

where the volume term depending on the function Cs is not conserved for any non-isentropic 
equation of state (12. ip . 

In even dimensions n = 2m > 2, 

[ pf{w/p)d^x (3.18) 
Jv{t) 

defines the enstrophy integral in terms of an arbitrary nonlinear function f of w/p, with 

= e • ((V A u) A ■ • ■ A (V A u)) = *(a;'") (3.19) 

" V ' 

n/2 times 

denoting the vorticity scalar of the fluid. (In the case when / is a linear function of w/p, the 
integral (I3.18P reduces to the trivially conserved circulation integral, which will be discussed 
further in the next section.) The transport equation obeyed by w is given by 

dtw = -V ■ (rou) - ma -w, w = *{u:"''^) (3.20) 

as obtained from (I3.19P and (13. 2p . To now evaluate the time derivative of the enstrophy 
integral (13.190 . we start by combining the transport equation (I3.20p and the Euler equation 
dH to get 

dtUJ = — u ■ Vw — mp^^w ■ cr, w = tjj / p, (3-21) 

and hence 

dt{pf{w)) + V ■ (up/(^)) = mf\w)w ■ CT. (3.22) 

Thus, in the case of isentropic compressible fluid flow, conservation of enstrophy holds due 
to (T = 0. In the general case of non-isentropic compressible fluid flow, we can rewrite the 
right hand side in equation (I3.22p via 

w ■ <T = -V ■ ((w ■ V^)e5) (3.23) 
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which arises from the curl expression (13. lip for the antisymmetric tensor cr and from the 
divergence identity 



V • w = *(V A LJ""-^) = (m - 1)(V A Lj) ■ *(w'"-2) = 0. (3.24) 
We thereby obtain, after integration by parts, 

VtT + V- (uT) + • ^ = -mesV ■ (w ■ VSf'{w)) (3.25) 

with 

T = pf{w), ^ = -mesf{w)wVS, (3.26) 

where / is a nonhnear function of w. Now, we find Ep^esV ■ {w - V S f (w)) = p^^ f"{uJ))Psw■ 
(Vzu A VS) ^ due to Ps ^ 0, and hence the right hand side in equation fl3.25p fails to 
be a spatial divergence since it is not annihilated by one of the spatial Euler operators. 
Consequently, for fluid flow with a non-isentropic equation of state, the enstrophy integral 
fl3.18p on any spatial domain V{t) transported in the fluid obeys 



^ / pf{w/p)d^x = -m [ esnw/p)w ■ (V5 A n)d" 

Jv(t) JdV(t) 

-ml esV ■ (w • VSf{w/p))(Fx (3.27) 
Jv{t) 

where the volume term that depends on the function is not conserved. 

Essentially, the underlying reason for non-conservation of both the enstrophy integral 
(I3.18P and the helicity integral (13. 4p for all non-isentropic equations of state is that es = 
J p~'^Psdp is necessarily a non-constant function of p. 

We now address the open questions of whether any generalizations of the helicity (13.40 
and enstrophy (13.180 are conserved for a non-isentropic equation of state or, more gener- 
ally, whether the non-isentropic fluid equations (I2.1l) - (l2.4p admit any conservation laws of 
vorticity form 

T(p,^,u,*(a;™)), m>l (3.28) 

where 

*K»)=(- = (3.29) 
1-07 n n = 2m + 1 

is the vorticity scalar in even dimensions and the vorticity vector in odd dimensions, respec- 
tively. 

Theorem 3.1. For any non-isentropic equation of state Ii2. 1\) the only nontrivial vorticity 
conservation laws jiS. 28\) for compressible fluid flow in dimensions n > 1 are given by 



d 
di 



[ f{S)wd''x= [ e5W ■ (VF(5) A n)t;"-V (3.30) 

Jv{t) JdV(t) 



for even dimensions n = 2m, where f{S) = F'{S) is a non-constant function of the entropy 
S . In particular, there are no special equations of state for which additional vorticity conser- 
vation laws are admitted in any even dimension, and no non-trivial vorticity conservation 
laws are admitted in any odd dimension n > 1. 
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The proof of this classification theorem will be given in Appendix |Bl 
The even-dimensional vorticity conservation law fl3.30p can be rewritten through integra- 
tion by parts via the identities 

CC7 = w = w ■ (V A u), V-w = 0. (3.31) 

This yields the equivalent conserved integral 

d 



dt 



[ w (u A Vf{S))d^x = [ {l\u\^Vf{S) - esVF{S)) ■ (w ■ n)rf"-V 

Jv(t) JdV(t) ^ 



(3.32) 



'V(t) JdVit) 

where n = 2m. We then see that the case f{S) = S has the form of a circulatory entropy 

^ / w ■ (u A V^)rf"x = / (J|u|2-e5^)w- (V^An)c/"-V (3.33) 
dt Jv(t) Jdv{t) 2 

which vanishes whenever the fluid flow is either irrotational, i.e. a; = and thus w = 
*{uj^~^) = 0, or isentropic, i.e. S =const. and thus VS = 0. Hence, in the general case, the 
conserved integral (13.321) can be understood to arise from conservation of circulatory entropy 
(I3.33P by the freedom to functionally redefine the entropy S* to S' = f{S). 

To conclude, we point out that the respective transport equations for the vorticity (13.291) 
in even and odd dimensions each yield locally trivial conservation laws (therefore falling 
outside of Theorem 13. ip in which 

T = *(a;™) = V ■ *(u Acu^-i) = dive (3.34) 

is a spatial divergence. In particular, in even dimensions n = 2m, we see the transport 
equation (I3.20p for zu takes the form of a conservation law (12. 8p given by the conserved 
density T = tu = V ■ (w ■ u), with the moving spatial flux ^ = — me^w ■ VS obtained 
from the identity (13.230 . We note this is precisely the local conservation law that underlies 
the conserved integral (I3.30p when the fluid flow is isentropic. Similarly, in odd dimensions 
n = 2m + 1 {m > 1), by combining the transport equation (13. 6p for zcr with the curl 
expression (13. lip , we obtain a conservation law (12. 8p where T = zu = V ■ (W ■ u) is a 
conserved vector-density and ^ = — u ® tu — me5W ■ VS is the moving spatial tensor-flux. 

4. Constants of motion 

A fluid conservation law in integral form (12. 9p on a domain V{t) C M" moving with the 
fluid yields a constant of motion 

^ I Td''x = (4.1) 
dt Jv{t) 

if the net flux through the moving domain boundary dV{t) is zero for all solutions of the 
Eulerian fluid equations ( ]2.2p -f l274l) . This condition has an equivalent formulation as a local 
continuity equation in the form (12.80 such that the moving flux vector ^ is divergence free 
in the domain V{t). Thus, a conserved density T determines a constant of motion (14. ip if 
and only if 

VtT + ■ (uT) = 0. (4.2) 

The resulting conserved integral (14.10 will be called a moving-domain constant of motion 
provided the conserved density is non-trivial, i.e. if T is not a spatial divergence div0 where 
the vector is a local function of t, x, u, p, S", and their x-derivatives. 

10 



Our classification of non-trivial kinematic and vorticity conservation laws in Theorems 12.11 
and 13.11 provides an immediate corresponding classification of moving-domain constants of 
motion of kinematic form fl2.24p and vorticity form fl3.28p - fl3.29p . 

Proposition 4.1. For compressible fluid flow in n > 1 dimensions, the only kinematic 
constants of motion 

^ [ T(t,x,u,p,S)rf"x = (4.3) 
dt Jv{t) 

on an arbitrary moving domain V{t) transported by the fluid consist of mass 



pd^'x (4.4) 

'Vit) 

and generalized entropy 

[ pf{S)d'^x (4.5) 
Jv{t) 

where f is any non- constant function of S. These constants of motion hold for all equations 
of state Hjy. 

In particular, no additional kinematic constants of motion arise for any special (non- 
isentropic or isentropic) equation of state or in any special dimension. 

Proposition 4.2. There are no vorticity constants of motion 
d f d f 

— I T{p, S,u,zu)d'^x = if n = 2m or — T(p, S,n,Tjj)d^x = if n = 2m + 1 

dt Jv{t) dt Jy(^i-^ 

(4.6) 

for compressible fluid flow with any non-isentropic equation of state Ii2.1\) in n > 1 dimen- 
sions. In the case of isentropic equations of state l{2.5\) . where 

e = j p-^P{p)dp {with es = 0) (4.7) 

is the internal (thermodynamic) energy density, the even-dimensional enstrophy 

pf{w/p)d''x {n = 2m) (4.8) 



given by any nonlinear function f yields the only vorticity constant of motion on an arbitrary 
moving domain V{t) transported by the fluid in n > 1 dimensions. 

In addition to the classifications given by these two Propositions, we can consider non- 
trivial constants of motion of the form 

^ / 0-nrf"-V = O (4.9) 

dt Jdv(t) 

holding formally for all solutions of the fluid equations f l2.2p -f E7^ . where dV{t) is the moving 
boundary of a domain V{t) transported by the fluid. Through the divergence theorem, any 
conserved integral (14. 9 p arises from a locally trivial conserved density T = div0 satisfying 
the moving-fiux condition f l4.2p in the domain V{t). In particular, a vector function yields 
a non-trivial conserved integral (14. 9 p if and only if 

A0 + u(V-0) = L'x-* (4.10) 
11 



holds for some antisymmetric tensor function ^ in terms of t, x, u, p, 5", and their x- 
derivatives. If is not identically divergence-free, the resulting non-trivial conserved integral 
(14. 9 p will be called a moving-boundary constant of motion. We remark that the corresponding 
local continuity equation fl4.1UI) can be viewed as a special type of lower-degree conservation 

law [iniE]. 

For isentropic fluid flow, a well-known example of a moving-boundary constant of motion 
is given by Kelvin's circulation theorem [I8j in n = 2 dimensions which states the circulation 
of the fluid velocity around any closed curve transported in the fluid is a constant of the 
fluid motion. The circulation integral has the form 

® u ■ ds = (p u ■ *nda = (p e-{uAn)da = — / zjd'^x (4-11) 
Jc{t) Jc{t) Jc{t) JV{t) 

related to the two-dimensional vorticity scalar vj = e ■ uj = *{W /\ w) (defined in terms of 
the spatial orientation tensor e in M^), where C{t) is any closed curve bounding a domain 
V{t) transported by the fluid in M^, *n is a unit tangent vector along C{t) or equivalently 
n is a unit normal vector, and da is the arclength element for C{t). As shown in [5], the 
circulation (14. lip has a generalization to all even dimensions n = 2m (m > 1): 

(f (uAa;'"-^)-dA= /" (u A o;™"^) ■ *nrf2'"- V = - /" wd'^'^x. (4.12) 

JdV{t) JdV(t) Jv{t) 

This integral defines a constant of the fluid motion for isentropic fluid flow, where dA = 
*hd'^'^~^a denotes the surface element for the moving-boundary hypersurface dV{t) in even 
dimensions n = 2m analogously to ds = *hda for moving-boundary curves in two dimensions. 

The natural question now arises as to whether the transport equations (I3.20p and (13.61) for 
the vorticity (13.291) in any even or odd dimension n > 1 yield a non-trivial constant of motion 
(14. 9 p when the fluid flow is non-isentropic. Our classification proof for vorticity conservation 
laws in Theorem 13.11 actually settles this question in a slightly more general form. 

Proposition 4.3. For compressible fluid flow with any non-isentropic equation of state Ii2. 1\} 
in n > 1 dimensions, there are no moving-boundary constants of motion Iji4-^ whose corre- 
sponding local continuity equation ^4-'^ for T = div0 is of vorticity form 28\) - f3. 29^) . In 
the case of isentropic equations of state \2.5\) . the only moving-boundary constant of motion 
with such a form is given by the generalized circulation Ili4.12\ ) in all even dimensions. 

In particular, for non-isentropic fluid flow, the vorticity transport equations written in the 
moving-boundary form (13.341) show that 



d 
di 



(f {u A Lj"^-^) ■ dA = m (f es{VSAiv"'~^)-dA (4.13) 

JdV(t) JdV(t) 



'dV{t) JdV{t) 

in even dimensions n = 2m, while in odd dimensions n = 2m -\- 1 
d 



dt 



{u A uj""-^) ■ dA = m (b es{VS Auj"'-^)-dA- (b u u"" ■ dA (4.14) 

dV{t) JdV{t) JdV(t) 



with dA = *n(i^"* ^cr, where ^ is a non-constant function of p for any non-isentropic 
equation of state (12. ip . 
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5. Hamiltonian correspondence 

A Hamiltonian formulation for the non-isentropic compressible Euler equations f l2.ip -fl2T 
in n > 1 dimensions (cf [ISj) is given by 

/ u \ / 5E/5VL \ ^ 

dA p \=n{ 5E/5p , E=ip|u|2 + pe(p,5)=p(i|u|2+ / p-2p(p,5)rfp), (5.1) 
\S j \ SE/SS J J 

in terms of the Hamiltonian operator 

/ p-i(V Au)- -V p-^VS \ 
%=[ -V- (5.2) 

V -{p-'^sy ) 

where E is the energy density of the fluid. This operator ( 15. 2p determines a Poisson bracket 

(5.3) 



. / 5G/5VL \ 

{jr^ g]^= h SF/Su SF/Sp SF/SS ) n 6G/6p 

\ 6G/6S J 



= j p-\\/ A u) ■ {5F/5u A 50 /5u) + 6G/6u ■ {V5F/5p - VSp-^5F/5S) 
- 6F/6u ■ {V6G/6p - VSp-^6G/6S) cTx 

satisfying (modulo divergence terms) antisymmetry and the Jacobi identity [T9], for arbitrary 
functionals T = j FdP'x and Q = J GdP'x where F and G are functions of t, x, u, p, S", and 
their x-derivatives. Here 5/5\i, S/Sp, S/SS denote variational derivatives, which respectively 
coincide with the spatial Euler operators E^, Ep, Es when acting on functions that do not 
contain time derivatives of u, p, S. 

Similarly to the isentropic case described in [5j, the Hamiltonian operator "H gives rise to 
an explicit mapping 

/ 6T/6n \ ^ f n\ /ryX 
-nl 6T/6p \ =X\ p \ = \ fj \ (5.4) 
V 6T/5S J \S J \4> J 

which produces infinitesimal symmetries X = t)J (9u + fjdp + (pds of the non-isentropic Euler 
equations from conserved densities T. In particular, these components of the symmetry 
generator 

f} = -p-\V Au)-5T/5u + VST/Sp-p-^VS5T/5S, f} = V-5T/5u, 4) = p-^VS ■ 5T/5u 

(5.5) 

will satisfy the infinitesimal invariance equations [191 



Dtfj + u • Vr) + r) ■ Vu - p'^fiVP + p-^V{Ppf] + Ps^p) = 0, 

(5.6) 

Dtfj + V ■ U17 + pV • r) = 0, Dt(f) + r) ■ + u • V</) = 

for all solutions of the non-isentropic Euler equations fl2.1l) -( |274l) . 

Evaluating the mapping fl5.4l) - fl5.5l) for the kinematic and vorticity conserved densities 
given by the conserved integrals in Theorems 12.11 and 13.11 we obtain the following two tables 
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of symmetries, all of which can be written as point transformations X = rdt + ^ J x + r/J 9u + 
'qdp + (pds given by 



f] = T] — rdfU — ^ ■ Vu, fj = 1] — rdtp — ^ ■ Vp, 



(5.7) 



where 77, rj, are functions of t, x, p, 5*, u, while r, ^ are functions only of t, x. 
For general non-isentropic equations of state fl2.ll) : 



Conserved Density T 


Conservation Law 


Symmetry X 


Description 


P 

pu 

pu A X 
p(tu — x) 
ip up + pe 
pS 

wS = *(a;™-i Au A V5) 


Mass 
Momentum 
Angular Momentum 
Galilean Momentum 

Energy 
Volumetric Entropy 
Circulatory Entropy {n = 2m) 




X A 9x + u A 9u 
t9x + du 
dt 




Nil 

Space translations 
Rotations 
Galilean boosts 
Time translation 
Nil 
Nil 


Here e is the internal (thermodynamic) energy density (\2.. 


7). Note the entropy (per unit 



mass) S can be replaced by an arbitrary function f{S). 
For polytropic equations of state fl2.25p : 



Conserved Density T 


Conservation Law 


Symmetry X 


Description 


tE — |p(u ■ x) 
i2^-tp(u-x) + ip|xp 


Similarity Energy 
Dilation Energy 


tdt + |xj 9x 
-^udu - \pdp 

— (tu — x)9u — ntpdp 


Similarity Scaling 
Galilean Dilation 


Here E = p(i up + e) = |p up + |/t(S')p"'^+^/" is the polytropic energy c 


ensity. 



As none of these symmetry generators X contain the pressure p, we see that the symmetry 
structure produced via the Hamiltonian mapping (15.51) for the non-isentropic Euler equations 
fi2.ip - fi2.4p is the same as in the isentropic case presented in [S]. This structure can be 
geometrically summarized as follows. 

Recall, the dilational Killing equation on M" is given by 

= r2g, Vl = const. (5.8) 

for a vector C,{^i where £ denotes the Lie derivative and g is the Euclidean metric tensor. 

Proposition 5.1. (i) For a general equation of state the Hamiltonian symmetries 

corresponding to the kinematic conserved densities for energy Ii2.15\) . momentum Ii2.12\) . 
angular momentum Ii2.13\) and Galilean momentum ^2.14^ comprise the generators of the 
Galilean group in n > 1 dimensions given by 



X = X 



n 



C\d^ + ntdt + -npdp 



(5.9) 



in terms of the solution of the Killing equation Ii5.8\) . (ii) The additional Hamiltonian sym- 
metries, consisting of the similarity scaling and Galilean dilation which correspond to the 
generalized energy densities Ii2.20\) and Ii2.21\) in the case of a polytropic equation of state 
Ii2.25\) . generate an extension of the Galilean group given by 



X = txJ5x 



1 ~ 77 ~ 

_nt% + -ntpdp 
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■1 A 

-tQu 

■2 



(5.10) 



where the vector x(x) is any irrotational solution of the dilational Killing equation = 
fig, VL = const., satisfying V A % = 0, on M". 

In contrast, the conserved densities for mass (12.110 . generalized entropy (I2.18p . and circu- 
latory entropy f l3.32p (or equivalently f l3.30p ) are mapped into the trivial symmetry X = 0. 
These conserved densities comprise all of the Hamiltonian Casimirs for the non-isentropic 
Euler equations ( I2.ip -( E^ in n > 1 dimensions. 

6. Gas Dynamics Conservation Laws 

The Euler equations fl2.ip -( p7^ for non-isentropic compressible fluid flow have an equiva- 
lent formulation in which the equation of state (12.1 p is inverted to give entropy as a function 
of pressure and density, S = S{p,p) with dS/dp ^ 0. The entropy transport equation (12. 4p 
then becomes 

dtp + u-Vp + F{p,p)V ■u = (6.1) 

in terms of the function 

(6.2) 

which replaces the equation of state (12. ip . For any specified choice of F(p,p), the entropy 
function S{p, p) can be recovered as a solution of the linear equation 

pSp + F{p,p)S^ = Q. (6.3) 

In particular, we can identify S'(j9, p) with the constant of integration arising from the char- 
acteristic ODE 

J-p = nv.p)lp. (6.4) 

Note that the simple state function 

p = 7p, 7 = const. (6.5) 

corresponds to a polytropic equation of state (12.251) which is equivalent to an entropy function 
of the form 5* = S{p/ p'^). 

The resulting formulation of compressible fluid flow thus consists of the dynamical equa- 
tions (16. ip . (12. 2p . (12.30 for pressure p, velocity u, and density p, respectively. These are 
precisely the equations of adiabatic gas dynamics [26j. Physically, the state function F{p,p) 
is then related to the sound speed c in the gas hj F/p = as seen from equation (16. 4p . 
The state function also determines the internal energy density e of the gas from the ther- 
modynamic relations de = —pd{l/p) = Cpdp + Cpdp combined with equation (16. 4p . which 
yields 

pep + F{p,p)ep = - (6.6) 
P 

for the function e{p, p). Note the homogeneous part of this linear equation (16. 6p involves the 
same linear differential operator pdp + F{p,p)dp appearing in the entropy equation (16.30 . so 
thus the general homogeneous solution of equation (16. 6p is given by f{S{p,p)) in terms of 
any particular solution of equation (16. 3p . where / is an arbitrary function of S{p,p). 

Since the dynamical pressure equation (16. 4p is related to the entropy transport equation 
(12. 4p by an invertible change of variables p = P{p, S) ^ S = S{p,p) (i.e. a point transfor- 
mation) when Sp ^ OT Ps ^ 0, there is consequently a one-to-one correspondence between 
the local conservation law structure of the gas dynamics equations (16. 10 - 06.20 . (12.20 . (12.30 
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and the non-isentropic Euler equations fl2.ip -( l2^ . as shown by general results in [HI [TO] . 
A complete classification of kinematic and vorticity conservation laws for gas dynamics in 
n > 1 dimensions now follows as a corollary to Theorems 12.11 and 13.11 

Theorem 6.1. (i) For a general state function F{p,p), all gas dynamics conserved densities 
of kinematic form T(i:, x, u, p, p) in any dimension n> 1 are given by a linear combination 
of the conservation laws for mass h2.11\) . momentum ^2.12\) . angular momentum l{2.13\) . 
Galilean momentum l[2.14 ), plus energy Ii2.15\) and volumetric entropy Ii2.19\) where e 



e{p, p) is any non-homogeneous solution of the thermodynamic energy equation lid. 6]) and 
S = S{p,p) is the general solution of the entropy equation 5^) . (ii) The only state function 
F{p, p) for which additional kinematic conserved densities arise is the polytropic case h6.^) 
with dimension- dependent coefficient'-) = 1 + 2/n. The admitted conservation laws consist of 
the similarity energy Ii2.20\) and the dilational energy \2.21\) given in terms of the polytropic 
energy density h2.23\) . 

Theorem 6.2. (i) For a general state function F{p, p), all gas dynamics conserved densities 
of vorticity form T{p,p, u, cc) = ifn = 2m orT{p,p, u, -aj) = ifn = 2m + l (with m > 1) 
consist of the conservation law for circulatory entropy Ii3. 33\) in even dimensions, as expressed 
in terms ofw = *{uj"'~^), VS = Sp(Vp — p~^F{p, p)'Vp), andcs = Cp/Sp, where S = S{p,p) 
is the general solution of the entropy equation Ii6. 3\) . e = e{p, p) is any non-homogeneous 
solution of the thermodynamic energy equation Ii6.6]) . and w is an antisymmetric tensor 
related to the vorticity scalar w by the identities VS. 31]) . (ii) There are no state functions 
F{p, p) for which additional vorticity conserved densities arise in any even or odd dimension 
n > I. 

7. Multipliers 

In index notation, the Euler equations for non-isentropic compressible fluid flow on M" in 
Cartesian coordinates [i = 1, . . . ,n) are given by 

+ J + p- V/ = 0, (7.1) 

Pt + ipu'),i = 0, (7.2) 

St + u'S,i = 0, (7.3) 

together with the equation of state 

p = P{p,S), Ps^O. (7.4) 

Similarly to the isentropic case studied in [5], any local conservation law T>tT +DiX^ = that 
holds for formally all solutions of these equations fl7.ip - fl7.4l) can be written in an equivalent 
characteristic form 

DtT + D,X' = {ul + u^u\, + p-'PpP: + p-'PsS:)Q, + {pt + {pu'),i)Q + {St + ,)i? (7.5) 
in terms of the multipliers 

Q, = E^,{T), Q = E,{T), R = Es{T). (7.6) 

In this formulation u'^,p,S are given by arbitrary functions of t and x*, with X* differing 
from X* by terms that are linear homogeneous in the fluid equations fl7.1l) - fl7.3p and their 
total spatial derivatives. Since the spatial Euler operators E^i, Ep, Es annihilate any spatial 
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divergence -Dj0\ there is then a one-to-one correspondence between nontrivial conserved 
densities T (modulo spatial divergences) and non-zero multipliers {Q\Q,R}. 

Necessary and sufficient equations for determining multipliers f l7.6p are given by apply- 
ing variational derivative operators 6/6u\ S/6p, 6/6S to the characteristic equation (17. Sp . 
yielding a linear homogeneous polynomial system in ul, pt, St, ulj, ptj, Stj, etc. whose co- 
efficients must separately vanish. The resulting determining equations for Q*, Q, R consist 
of the adjoint of the determining equations for symmetries (15. 6p . plus additional integrablity 
equations (Helmholtz conditions) for Q\ Q, R to have the variational form (17. 6p . Thus, 
multipliers can be characterized as adjoint-symmetries that have a variational form |1| [TO]. 
Moreover, the corresponding conserved density T (up to an arbitrary spatial divergence) 
can be constructed explicitly from the multipliers Q*, Q, R by means of homotopy inte- 
gral formulas fL9[ HJ [10] or by an algebraic scaling formula [1] based on invariance of the 
fluid equations under dilations t — )■ At, — )■ Ax*. The determination of multipliers and 
hence of conservation laws is thereby reduced to an adjoint version of the determination of 
symmetries. 

To conclude, in the following two tables we list the multipliers for, ffistly, the non-isentropic 
kinematic conservation laws (I2.1ip - (l2.15p and the non-isentropic vorticity conservation laws 
(I3.30p and (equivalently) (I3.32p -f l3.33p . all of which hold for a general equation of state (12. ip : 
and, secondly, the extra non-isentropic kinematic conservation laws (I2.20p - (l2.2ip holding 
only for the distinguished polytropic equation of state (12. 25 p . 



Conserved density T 


Description 


Qi = 5T/5u' 


Q = ST/6p 


R = 6T/6S 


P 


Mass 


1 








h 

pu 


Momentum 


P^^ 







p{u^x^ — u^x^) 


Angular momentum 


p{x^5i - x-'^f ) 


u^x^ — u^x^ 





p{tu^ — x^) 
pilu^Uk + e) 
pS 
zuS 


Galilean momentum 


pt5^ 


tu^ — x^ 





Energy 
Volumetric entropy 
Circulatory 


pUi 




\u^Uk + e + pCp 
S 



pes 
P 




entropy {n = 2m) 









Conserved density T 


Description 


Qi = 6T/6u' 


Q = 6T/6p 


R = 6T/6S 


tE — ^pu^Xk 
t^E - p{tu^ - \x^)xk 


Similarity Energy 
Dilational Energy 


pitui - ixi) 

tpitUi - Xi) 


^u'^ituk - Xk) 

+t(l + f )/€(5)p2/- 
\{tU^ - X^){tUk - Xk) 

+t\l + ^)k{S)p"^ 


tes 

ees 



Here E = pi^u^Uk + e) = \pu^Uk + |k(S')p^+^/" is the polytropic energy density, in terms 
of the internal energy density e = ^k(S')p^/". 



Appendix A. Proof of Theorem 12.11 and Proposition 14.11 

The time derivative of a kinematic conserved density r(t, x*,p, m*, S) is given by 

AT = -Tsu'S^, - Tpiu'p,, + pu\i) - T,.{u^u\, + p-'P;) + Z. (A.l) 

By applying each Euler operator Ep, Es, E^i to ( lA.ip . we get three linear inhomogeneous 
expressions in u\j, u\i, p,\ S*/, whose coefficients must separately vanish. This yields the 
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system of determining equations 

p-'PsTu^,, + Ts^^u' + Tts = 0, (A.2) 

pTpxi + T^i^jU'' + Ttui = 0, 

T,,.n* + p-^PpT^^x^ + Ttp = 0, (A.4) 

{pTps - Ts)5,j - P^'PsTu^u. = 0, (A.5) 

P Tpp5ij — PpT^iy^j = 0, (^-6) 

{pT„,p - T^.)Ps - pPpT^.s = 0, (A.7) 

pTu'p - T^^ = 0, (A.8) 

to be solved for T. Here 6ij = 5(^ij) ■v^ g denotes the Cartesian components of the Euchdean 
metric tensor on M". 

We start by integrating equation flA.Sp with respect to p and u*, yielding 

T = pf{t,x\u\S)+g{t,x\p,S). (A.9) 

Now substituting (lA.Qp into flA.7p we get fu^s = ^ and then integrating gives 

f = fit,x\u') (A.IO) 

where an integration constant c(t,x\S) has been dropped since we can absorb pc(t,x\S) 
into the term g(t, x\ p, S) in (lA.Op . 

Next, substituting (IA.9P and (lA.lOp into both (IA.6D and OA. 51) . we obtain 

pgppSij - UuiPp = 0, ( A. 1 1 ) 

{P9ps - gs)^ri - = 0. (A. 12) 

Separation of flA.lip with respect to yields the equations 

fu^ui = c{t,x^)5ij, (A. 13) 

gpp = c{t,x'')p-^Pp, (A.14) 

where c{t^x^) is a separation constant. Substituting (1A.13P into (lA.12p . we get 

pgps-9s = c{t,x'')Ps. (A. 15) 

By writing flA.121) and flA.151) as {pgp — g — cP)p = {pgp — g — cP)s = 0, then integrating 
with respect to p and S, we obtain a first-order linear equation with respect to p, 

pgp-g = cP-c,{t,x'). (A.16) 

Integration of flA.16p and flA.lSp then yields 

f = Cj{t,x')u^ + lc{t,x')u^Uj, g = p~^ci{t,x') + c{t,x') j p''^P{p, S)dp, (A.17) 

after integration constants have been absorbed as before. Consequently, from (1A.9P we have 

T = ci(t, x') + C2(t, x\ S)p + Cj{t, x')pu^ + c(t, x'){\pu^Uj + pe(p, S)) (A.18) 

which gives the general solution of the determining equations flA.8p . (lA.7p . (lA.6p . flA.5p . with 
e(p, S) given by (12.170 in terms of P(p, S). Note the term ci in (lA.lSp is a trivial conserved 
density, and so we will put ci = 0. 
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By substituting (lA.lSP into the remaining determining equations (lA.4p . (lA.Sp . (lA.2p . we 

get a system of equations to be solved for c(t, x*), Cj(t,x^), C2(t,x\ S), ci(t,x^): 

c^i = 0, (A.19) 

(A.20) 
(A.21) 
(A.22) 
(A.23) 
(A.24) 



Cjxi + Cja-i + Ct^ij — 0, 
C2xi -\- Cit =0, 
Cix.Pp + Cte + CtpCp + C2t = 0, 

Cix.Ps + ctpes + pC2st = 0, 



To begin, we note flA.lQp implies c = c{t), and hence flA.20p has the form of a time- 
dependent dilational Killing vector equation on Cj. Its solution consists of a linear polynomial 
in given by 

c^ = Co^{t) + Cuj{t)x^ ' \d{t)x,, Cujit) = -Ci,,{t). (A.25) 

Then, using the trace of flA.20p in (IA.22p . followed by differentiating with respect to x*, we 
find 

C2tx^ = 0. (A.26) 

(A.27) 



Integration of flA.24p and flA.26P gives 

C2=C2{t,S) + C2{x'). 

Substituting ([OTI) and ([OSl) into f[OTl) . we get 



Cq^^ C-\j^X'^ ~t~ Q c x^. 



^lij-^ ~r 2^ "^i- (A. 28) 

By differentiating with respect to x^ and antisymmetrizing in i,j, we find C[^j = 0, so thus 
Ciij is constant. Integration of ( 1A.28I) thus yields 

(A.29) 



C2 



C2 C'gjX -|- XiX 



where C2 is a constant. Then, differentiating ( 1A.29P with respect to t, we get — Cq^x* + 
ic'"xjX* = which splits with respect to into Cq^ = 0, c'" = 0. Hence we obtain 

Coi = CLoi + diit, c = bo + bit + 62^^, (A. 30) 

and, from fCOSjl and (1X29]) . 

Ci = aoi + aiit + CiijX^ - ^{bi + b2t)xi, C2 = C2 - aux' + ^b2x'xi, (A. 31) 

where aoi,aii,bo, bi, 62 are constants. Now, through (]A.3ip . we can rewrite (]A.22p and ( IA.23P 



as 



inP + pe)^ 



-(pc: 



2t)p, 



\nP + pe)<. 



■(pc: 



'2t)S- 



(A.32) 



Integrating flA.32p with respect to p and S*, we obtain 

(pe - inP)c' + pc2t = a{t). (A.33) 

By considering the separation of (1A.33P with respect to p, we have the following two cases. 

Case c' = 0: Hence c = const. = 60 and 61 = 62 = in (lA.30p . Now splitting (]A.33P with 
respect to p yields a = and C2t = 0, so thus 

52 = b{S) (A.34) 
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where, through (lA.27p . we have absorbed an integration constant into C2. Hence, ( lA.lSP 
reduces to 

T = C2P + b{S)p + aoipu' + aiip{tu' - x') + Cujpu'x^ + boHpu'ui + pe{p, S)) (A.35) 

with arbitrary constants Ooi, a^, Bq, Cuj = —Ciji, C2, and an arbitrary function b{S), where 
e{p,S) is given by fl^TT]) . 

Case c' 7^ 0: Dividing (1A.33P by p, then differentiating with respect to p and multiplying 
by p^, we get p^c'^Cp — ^n{p~^ P) p) = —a which separates with respect to t and S into the 
equations 

p\ep - ln{p-'P)p) = d, (A.36) 
c'd = -a, (A.37) 

where d is a. constant of separation. Substituting e from (12.1 7p into (lA.36p . we obtain a first- 
order hnear equation with respect to p, —^npPp + (1 + ^n)P = d, whose general solution is 
given by 

P{p, S) = K(^)p^+2/n ^ ^/^^ ^ ^/2). (A.38) 

Since we can (without loss of generality) shift the pressure P by an arbitrary constant, we 
put d = 0. Then, flA.37p yields a = 0, while (I2.17P becomes 

eip,S) = lnp-^P + i{S). (A.39) 

Substituting flA.39p into flA.33p . we get c'e{S) + C2t = which gives 

cB + C2 = b{S). (A.40) 

As a result, (lA.lSP becomes 

T =C2p + b{S)p + aoipu' + aiip{tu' - x') + Cujpu'x^ + &o(|pm'm* + |nK(5)p^+2/") 

+ bi{t{\pu'u^ + in/t(S)pi+2/") - \px'ui) + b2{t^{\pu'u^ + in/t(S)pi+2/") - tpx'u^) 

(A.41) 

with arbitrary constants aoi, aii,bo, bi, 62, Cuj = — Ciji, C2, and an arbitrary function b{S). 
Also, (1A.38P and (1A.39P respectively reduce to 

P(p, S) = k{S)p'^'/\ e(p, S) = in«:(5)p2/". (A.42) 

Appendix B. Proof of Theorem 13.11 and Propositions 14.21 and 14.31 

To begin, we write out the component form of the fluid curl, the vorticity scalar and 
vector, along with their transport equations: 

oj ^ = \{u\^ - u^^) = u^\^\ (B.l) 

^ = e*i^'i-*™^™Wi,ji ■ ■ ■ Ui^j^, m = n/2, (B.2) 

zu^w' = e'^^'^^-^-'-^tUj.k, ■ ■ ■ coj^k^, m = {n- l)/2, (B.3) 

= {ukuj^')^^ - {ukuj^'): = Uk:uj^'' - Uk^u^"" - wV^fc, (B.4) 

vjt = -{vju^ + p^^w''^Pj)^i = -vju\i - u'-w^i + p^'^Psw^^p^jS^k, (B.5) 
w\ = {w^u' - w'v? - p-^W'^''P^k),j = ^^u\j - zo'u^^j - v?w\^ + p-'^PsW'^^p^jS^^. (B.6) 
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In addition, we will need the component form for 

w o w'^ = e^^-^^^--*--^--a;.,,, ■ ■ ■ m = n/2, (B.7) 

W o W^^' = e'^^^-'--^--'--uj,^, . ■ m = (n - l)/2, (B.8) 

VU7 O = W^''ujjk,i, VtU ^ G7-'',i = W^^^UJkl,i, (B.9) 

along with the identities 

= 0, w'{j = 0, nui^^Wfcj- = zu5i (B.IO) 

lyiO'fc) = |yfe> = 0, PF^^'^fc = 0, (B.ll) 

w\i = 0. (B.12) 

Here e*^"'*" = e'*^"'*"' f-)- e, = ^ g are the Cartesian components of the spatial orien- 
tation tensor and the Euclidean metric tensor on M"; round brackets denote symmetrization 
of the enclosed indices, and square brackets denote antisymmetrization. 

We proceed by explicitly solving the determining equations (12.101) for conserved densities 
of vorticity type in even and odd dimensions n > 1. 

Case n = 2m: For a conserved density of the form T{S, pjU"^ the time derivative is 
given by 

VtT = D,{-u'T) + u\,A - p-\PpP: + PsS:)T^^ + p-^Psw'' p^^S^^T^ (B.13) 
after use of identities fIB.lOl) . where 

A = T - pTp-wT^. (B.14) 
To begin we substitute flB.131) into the determining equations 

= Es{VtT) =u'^p'^PsipT^^u^ + mpl'wikT^j^) + p' p''^Ps{pT^^p + |a7T„»„ - - fT^^s) 
+ w^p~'^Ps{pTu^.^ - mwijp^^T^^) + tiuAs (B.15) 

and 

= E,{VtT) =u'^p'^PsUTu^u^ - mS'w^kT^,^) + S^p-^Ps^T^^ - \wT^.^ - pT^.^) 

+ zu'p~'^PsifT^^^ - mWijS^^T^^) + truAp (B.16) 
where we have introduced the notation 

tiu = 6iju'^ = u\i, u'^ = l{u\^ + u{') = / = pPp/Ps. (B.17) 

Since T does not contain any derivatives of w, S and p, the coefficients of the separate terms 
w^S^^ , ■^/p/) vj^ in both flB.151) and flB.16p must vanish, yielding 

Tt^vj = 0, T^ijjj = 0. (B.18) 

Integrating fIB.lSp . and dropping a kinematic term c{S,p,u^) that does not involve zu, we 
get 

T = a{S,p)w. (B.19) 

Then f lB.lSP and f IB.lGp reduce to = tiuAs = tiuAp, and hence we have = As = Ap 
which gives 

ipap)s = {pap)p = 0. (B.20) 
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Integration of these first-order linear equations for Op yields 

a = b{S)+c\np, c = const.. (B.21) 
Thus from flRTOjl and I^B2T\f we have 

T = b{S)rz! + cw\np (B.22) 

and 

VtT = Di{-u'T) - cwu\i + {h{S) + c\Tip)p-'^Psw'^p^iSj. (B.23) 
The third term in flB.23p can be written as a spatial divergence —Di{B{S,p)w^^pj) = 
w^^ Bsp^iSj through the identities fIB.lOp . where B{S,p) satisfies Bs = p^^{b{S) + c\np)Ps. 
Thus, flB.23p becomes 

VtT = -Di{u'T + Bw'^p^j) - cwu\i. (B.24) 
Now we substitute flB.24p into the final determining equation 

= Eui(VtT) = czu^i + mctru^^Wji (B.25) 

Since w has no dependence on u"^^ (which is linearly independent of a;*-'), we obtain 

c = 0. (B.26) 

Therefore, flB.26P and flB.22p give the result 

T = b{S)zu (B.27) 

with 

VtT = -Di{u%S)w + B{S, p)w'^p,j), B{S, p) = p-^ J b{S)PsiS, p)dS. (B.28) 

We now note that if 6 = const, then flB.27p is a trivial conserved density 

T = bw = b{w'^Uj)^i = Di{bw'^Uj). (B.29) 

Its corresponding moving- flux ^* = bp^'^Pw^^ pj is non- vanishing due to the form of (]B.28p . 
This completes the proof in the even- dimensional case. 

Case n = 2m + 1: For a conserved density T{S, p,u^,uj^), similarly to the previous case, 
its time derivative is given by 

VtT = D,{-u'T) + u\,{T - pTp - zo'T^,) - p-^PpP,' + PsS:)T^u^ + p-^PsW''''p.,S,kT^.. 

(B.30) 

The determining equations = EsiVtT) and = EpiVtT) then yield 
=w\^p-^Ps{pT^.^^ + p/W/T„.^.) + p-^PsP:{pT^^p - T„,: - fT^.s - VT., 

+ u'^w^T^^s + p-'PspTu^u^ + p-'Psp!Wii''T^.^,) + tiu{Ts - pT^s - w'T^.s) (B.31) 

and 

=w\^p-^PsUTu.^^ - S,'W/T^.^,) + p-^PsS^fTu^s + - pT^p + W^., 

+ u^'izu.T^.p + p'^PpT^.u^ - p-^PsS^Wa^'T^^u^) - truipT^, + zu%^,) (B.32) 

using the notation (lB.17p . Since T does not contain any derivatives of w\ the first term in 
both (lB.3ip and (1B.32P must vanish modulo the identity (lB.12p . This implies 

^roiro* = 0, (B.33) 
T^,^, =aip,u^ui'')5ij. (B.34) 
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Applying the derivative operator (9„fe to (1B.34|) . antisymmetrizing in [jk], and taking the 
trace over (ij), we get (n — l)a„fe = 0. Hence, in n > 1 dimensions, 

a^k = 0. (B.35) 

By also applying the derivative operator d^k to (lB.34p . we similarly get 

a^k = 0. (B.36) 

Integration of flRMj) . (EM]), flRM]) . (^36^ then yields 

T = a(5,p)u*o7i + 6i(5,p)ii7\ (B.37) 

Here we have dropped an integration constant c{S,p,u^) since it does not involve -gj* (i.e. it 
is of kinematic form). 

The determining equations (lB.3ip and (1B.32P thereby reduce to 

= u'^BisWj - ptTuw'BiSp + p^Wip-^A, (B.38) 

= u'^BipW, - tiuw\pBipp + B,p) - S^Wip-^A, (B.39) 

with coefficients 

A = {pap - fas - (m + l)a)Ps, (B.40) 

Bi = aui + hi. (B.41) 

In (]B.38P and (IB. 390 the respective coefficients of the terms and must vanish, which 
yields A = Since P5 ^ 0, we get 

= pap - fas - (m + l)a. (B.42) 

Then, since the coefficient of u^^ must vanish in both (1B.38P and (IB. 390 . we obtain 

^{jBi)s = w^pBkps^ij-, (B.43) 
'co(^jBi)p = zu\pBkp)p6ij. (B.44) 

By taking the product of (1B.43P with wiWh and antisymmetrizing in [ih] and [j7], we find 
pw^ Bkps'OJ[h5i][j'(Ui] = 0. This implies Bkps = 0. The same antisymmetric product applied to 
(]B.44p then implies {pBkp)p = 0. Hence (IB. 430 and (]B.42p become zu(^jBi)s = Wf^jBi-^p = 0, 
which yields 

B^s = B,p = 0. (B.45) 
From (IB. 411) and (IB. 421) we thus obtain 

a = 0, b' = const.. (B.46) 

Hence (IB.37P reduces to 

T = biw' = Dj{biW'^''uk) (B.47) 
which is a trivial conserved density, with 

VtT = D,{b,{2u^'w'^ - W\p-^),k)) (B.48) 

as given by (]B.30P combined with the identities u\jzu^ = Djiu'w^) and p^'^PsW^^^ p^S^k = 
-W'^^{p~^)jP^k = Dk{-W'^^P{p~^).j) obtained via (IRT2|) and (IBnTj) . The corresponding 
moving-flux has the form ^* = —bj(u^w'^ + p^'^W^^^ p^k) which is non- vanishing. 
This completes the proof in the odd-dimensional case. 
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